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Randomly diluted quantum boson and spin models in two dimensions combine the physics of 
classical percolation with the well-known dimensionality dependence of ordering in quantum lat- 
tice models. This combination is rather subtle for models that order in two dimensions but have 
no true order in one dimension, as the percolation cluster near threshold is a fractal of dimen- 
sion between 1 and 2: two experimentally relevant examples are the 0(2) quantum rotor and the 
Heisenberg antiferromagnet. We study two analytic descriptions of the 0(2) quantum rotor near 
the percolation threshold. First a spin-wave expansion is shown to predict long-ranged order, but 
there are statistically rare points on the cluster that violate the standard assumptions of spin-wave 
theory. A real-space renormalization group (RSRG) approach is then used to understand how these 
rare points modify ordering of the 0(2) rotor. A new class of fixed points of the RSRG equations 
for disordered ID bosons is identified and shown to support the existence of long-range order on 
the percolation backbone in two dimensions. These results are relevant to experiments on bosons 
in optical lattices and superconducting arrays, and also (qualitatively) for the diluted Heisenberg 
antiferromagnet La2 ( Zn ,Mg) a, Cui _ ^ O4 . 

PACS numbers: 05.30. Jp,64.60.Ak, 74.81.-g,75.10.Jm, 75.10.Nr, 75.40.Cx, 75.40.Mg 



I. INTRODUCTION 

The physics of classical lattice percolation has been 
used to model a great variety of physical systems with 
considerable success and recent results on 2D percola- 
tion have put this classical theory on a firm mathematical 
footing^ Several materials of current interest are well de- 
scribed by combining random dilution with lattice quan- 
tum degrees of freedom, as in a magnetic material in 
which some quantum spins have been removed by chem- 
ical dilution. 3 Aside from its interest as a microscopic 
description of dilution in materials, percolation is impor- 
tant as the simplest nondeterministic process for gener- 
ating a "fractal" , an object of fractional dimensionality: 4 
as explained later in this introduction, there is a geomet- 
ric phase transition in randomly diluted lattice systems 
with a fractal structure at the transition point. 

The focus of this paper is on models for randomly di- 
luted interacting bosons. These models describe super- 
conducting Josephson-junction arrays with only onsite 
interactions and a finite density of defective junctions, 
or bosons in optical lattices where some sites remain 
empty. Percolation of superconductivity at short length 
scales has also been discussed as a model for inhomogene- 
ity observed by scanning tunneling microscopy on the 
high-temperature superconductor BSCCO. 5 Finally, the 
0(2) rotor is qualitatively similar to the nearest-neighbor 
s = 1/2 antiferromagnetic Heisenberg model in that both 
have long-range order in two dimensions but only alge- 
braic correlations in one dimension. The diluted Heisen- 
berg model has been studied numericall y 316 : 7 and in neu- 
tron scattering experiments on La2(Zn,Mg) x Cui_ a: 04 3 . 

Our interest in interacting bosons suggests that we con- 



sider the diluted 0(2) quantum rotor model. More gen- 
erally, the O(N) quantum rotor model^ has an A com- 
ponent unit vector hi living at every lattice site i. The 
effect of dilution is incorporated by defining the symbol 
ctij to be unity if the bond ij is present in the diluted 
system, and zero otherwise. For concreteness, we con- 
sider bond dilution though the results are more general. 
Then, the Hamiltonian for a particular realization of the 
dilution {cxij} is given by: 

H =^Y^ L 1 - J Yl "'J" 1 ' h o (!) 

i <ij> 

where, U and J are positive coupling constants, and the 
components of the angular momenta \L] a p are given by 
[L] Q /3 = Panp - ppn a (a ^ [3) where [n a ,pp] = iS a p, 
and the greek indices take values from 1 to A. Thus 
the kinetic energy part of the Hamiltonian is just L 2 = 
J2 a <(3 Lap- For undiluted lattices of dimension d > 2, 
the quantum O(N) rotor model has long-range order, 
as long as quantum fluctuations (measured by the ra- 
tio U/J) are below some non-zero, A-dependent, critical 
value (U I J) c . However, a chain of quantum rotors has no 
long-range order for any A > 2 unless U = 0. 8 

In a certain sense, a planar material with a large frac- 
tion of diluted bonds (or sites) has geometry between 
one- and two-dimensions. When the dilution fraction p 
reaches a lattice-dependent critical value (p c — 1/2 for 
square lattice bond dilution), there is a geometric phase 
transition: for p < p c there is an infinite nearest-neighbor 
connected cluster of undiluted sites for essentially all re- 
alizations of dilution, while for p > p c there is no infinite 
cluster i At threshold (p = p c ), the number of bonds of 
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the cluster, N(r), contained in a small circle of radius r 
around a given bond of the cluster: N(r) ~ r d f , with the 
"fractal" dimension d f = 91/48 « 1.896[Ref. Q- Now, 
an undiluted chain of sites has N(r) ~ r and an undiluted 
square lattice, iV(r) ~ r 2 , so, in this sense, the cluster 
is "between" one- and two-dimensions. Note, however, 
that other geometric properties of the cluster have other 
dimensionalities, some, in fact, less than one, 1 

Allowing quantum fluctuations (U ^ 0), raises an im- 
portant question about the nature of order in diluted, 
planar systems: Does the order on the critical cluster 
which exists for U — persist for non-zero U? Early 
numerical and analytic work on the diluted Heisenberg 
model suggested that long-range order vanishes before 
the percolation threshold;*!*^ although recent numerics 
disagree£ii This paper begins to resolve this question 
analytically by coupling a spin-wave approach to a com- 
plementary real-space renormalization group (RSRG) 
analysis. 

The spin-wave approximation, described in Section II, 
provides a natural first step for describing the effects of 
weak quantum fluctuations on the O(N) rotor model at 
threshold. We show that the fracton dimension d s of 
the cluster, which is close to 4/3 for percolation in any 
dimension^ is the relevant dimension to consider when 
discussing ordering near percolation threshold. Since the 
fracton dimension is greater than one, quantum fluctu- 
ations of the order parameter on an average site of the 
cluster are small, according to our computation. We also 
cite a rigorous result that implies, within the framework 
of the spin-wave approximation, that all sites have small 
quantum fluctuations provided that p < p c . 

In principle, this is not enough to argue that the super- 
fluid order is stable to quantum fluctuations. One might 
worry whether the stability to quantum fluctuations per- 
sists beyond the weak quantum fluctuations allowed by 
the spin-wave approximation. Further, at p = p c where 
the spin-wave approach only describes the behavior of 
an average site, one might be concerned that quantum 
fluctuations suppress superfluid order on a set of sites of 
measure zero on the percolation cluster. If these special 
sites are such that in their absence the full connectivity 
of the cluster is lost, then there still may be no true long 
range order. For instance the order may be destroyed on 
the long ID segments connecting the large blobs on the 
percolation cluster's backbone (Fig. QJ. 

In Section III, we identify an effective model for ro- 
tors on the cluster backbone that incorporates fluctua- 
tions along long links beyond the linear approximation. 
The physical idea, which we develop in detail is that 
these fluctuations occur on a qualitatively faster time- 
scale than fluctuations in the blobs. We refer to this as 
the slow blob approximation or SBA. Within the SBA, 
we compute a low-energy effective Hamiltonian for the 
O(N) rotor model on the percolation cluster. Due to the 
one-dimensional topology of the backbone, the effective 
Hamiltonian is a one-dimensional O(N) rotor model with 
strongly- inhomogeneous charging {ui} and exchange {ji} 



energies: 



Hsb = Y L ? ~ ■ ( 2 ) 

i i 



where the coupling constants {ji}, {ui} are independent, 
random variables, with distributions P(j) and R(u) re- 
spectively. Like the in lJT]l. these coupling constants 
are randomly distributed. However, the distributions 
P(j),R(u) are continuous for small j <C J and small 
u -C U, respectively. In contrast the ot^ have a discon- 
tinuous, bimodal distribution. Note that the SBA drasti- 
cally simplifies the geometry of the problem. In the SBA, 
the effective model © is truly one-dimensional: namely 
the ID O(N) rotor model with strong disorder. 

To reach this effective description, we focus on the 
backbone of the incipient infinite cluster. Most bonds lie 
off the backbone: our physical idea is that these "dan- 
gling" bonds will exhibit long-range order if and only if 
the backbone bonds do so. This is reasonable, since at 
[7 = 0, the dangling bonds play no role in communicat- 
ing phase correlations. Moreover, previous work follow- 
ing this line of analysis successfully explains a number 
of basic phenomena of ordering near percolation thresh- 
old. For example, the finite temperature phase diagram 
in any dimension of classical magnets diluted with non- 
magnetic impurities (e.g. Rb p Mni- p F2 which has a large 
5 = 5/2 local moment), follows from analyzing back- 
bone thermodynamicsii2ii£ Further, for quantum mag- 
nets with discrete symmetry (such as the transverse field 
Ising model) , backbone physics accounts for the zero tem- 
perature phase diagram as well^ The present work sim- 
ply extends the previous treatments to cover continuous, 
quantum degrees of freedom. 

Recently, the real-space renormalization group ap- 
proach to quantum systems with strong disorder, per- 
fected by D.S. Fisher in the context of random-exchange 
Heisenberg antiferromagnetic spin chains^ has been ex- 
tended to the ID 0(2) rotor model with strong disorder 
by E. Altman et.alii& In section IV and V, the machin- 
ery of this approach is brought to bear on the effective 
model in Eq|21 The result indicates that for the 0(2) 
rotor, superfluid order persists at threshold even with 
a large amount of quantum fluctuation, namely, up to 
U "/ J ' « {U/J)kt ~ tt 2 /4, where, (U / J)kt is the location 
of the Kosterlitz-Thouless transition of the clean ID 0(2) 
rotor model. Section V treats the range U/ J < {U/J)kt, 
while section IV discusses the loss of superfluid order 
which occurs just above (U/J)kt- 

Section VI concludes the paper with a discussion of 
the phase diagram of the 0(2) rotor model with quenched 
dilution. Two appendices present details of the real-space 
renormalization group computations. 
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II. SPIN WAVES ON THE CRITICAL 
PERCOLATION CLUSTER 



There are two main results of this analysis: First, ex- 
actly at threshold (p — p c ) and except for certain rare 
points ( "geometric fluctuations" ) , long-range order is sta- 
ble to weak quantum fluctuations. Second, away from 
threshold, on the percolating side (p < p c ), long-range 
order is stable to weak quantum fluctuations on any site, 
for almost all clusters. 

Clearly, if U — in we are in the classical limit of 
this model, and the ground state simply has hi = const. 
on all sites of a connected cluster. The question that 
we will address first is, exactly at the percolation thresh- 
old, does this long range order on the percolation cluster 
survive the addition of small quantum fluctuations? We 
will address this question within a spin wave calculation. 
Since the spin wave calculation is essentially identical for 
all N > 2, we specialize below to the 0(2) model: 



H = — 



2 



J ^ ati. 

<ij> 



COS 



(3) 



parameter is readily calculated: 



< expifa > 



exp(-- < fa 2 >) 



<fa > 



du> 



uj 2 1 

U 



JT 



(7) 
(8) 



Note that < faf > is formally of order y/U /J. For 
small U/J, it will therefore be small so long as [T - ^]^ 
is finite, and the long range order will survive quantum 
fluctuations. Below we will establish the required finitc- 
ness. 

Eqn. IHlmay be usefully rewritten in terms of the eigen- 
values and eigenvectors of the matrix T. Since T is a 
real, symmetric positive semi- definite (x T Tx > 0) ma- 
trix, it has real eigenvectors and eigenvalues (xn, A„) with 
A rl > (the zero eigenvalue is obtained for the uniform 
vector on connected sites): 



(9) 



The fa represent the phase of the bosons at site i. The 
operator rij has integer eigenvalues and physically repre- 
sents the excess boson number at each site. The fa and 
ni on a site are conjugate variables: 



[m, expifa] — Sij expi(/ 



(4) 



Clearly, this is just the eigenvalue problem for the Lapla- 
cian on the lattice defined by the a%j, which in this case 
is taken to be at the percolation threshold. We get 



„ V An 



(10) 



This model may also be taken to represent an array of 
Josephson junctions. U is a measure of the charging en- 
ergy that induces quantum fluctuations of the phase, J 
is the strength of the Josephson coupling. 

In the absence of quantum fluctuations (U = 0), the 
classical ground state is simply fa = const. The effect of 
turning on a small U/J may be addressed in a harmonic 
( "spin- wave" ) approximation by simply expanding the 
cosine: 



(•5) 



It is easier to work with the equivalent Euclidean 
(imaginary time) action: 



= 2 duj (~jj 6 v + JTij)fa(uj)fai(-u}) (6) 



where we have rewritten the potential energy part using 
Tij = c<ij(25ij — 1), and the frequency integral has dm — 

duj 

In this quadratic approximation, the superfluid order 



The eigenvectors are assumed to be normalized on the 
percolation cluster: 



E *n(j) 

j Gperc. cluster 



l 



(ii) 



To argue for the finiteness of the integral in Eqn. i|l(J[l , 
we first consider the site averaged value of the RMS phase 
fluctuation: 



< 



>Pe 



w r E <*<> 

y iG P 'ere. cluster 

N r ^ .r\ 



A,, 



(12) 
(13) 



Here N p is the total number of sites on the incipient 
cluster and we have used Eq. Q10[l. 1)11(1. In terms of the 
density of states per site p(X) for the eigenvalues of the 
operator T 



KA) = ^$>(A-A„) 



we have: 



< <p > Pe 



dX 



(14) 



(15) 
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Finiteness (or lack thereof) of the phase fluctuations is 
determined by the behaviour of the density of states at 
very small A. In the thermodynamic limit, the asymp- 
totic behaviour of p(X) is known!, since T is proportional 
to the Laplacian on the percolation cluster. We have 

p(A) ~ A^ -1 (16) 

where d s is called the fracton dimension. (On a regular 
one dimensional line, we note that d s = 1, which gives 
the well known logarithmic divergence in the RMS value 
of the phase fluctuations). For percolation clusters in 
spatial dimension d = 2, d s = 1.32 (sec Ref,i). In fact 
for percolation clusters in all dimensions, its value is very 
close to 4/3. 

Thus, for percolation in d — 2 we have: 

< 2 >Perc~ J d\\-° Si (17) 

which is clearly finite when integrated over small A. Since 
the fracton dimension of a percolation cluster does not 
vary much with spatial dimension, this finiteness of the 
site averaged RMS phase fluctuations holds for percola- 
tion in any dimension d > 2. 

To reinforce this point, we turn to an alternate repre- 
sentation of the suppression of superfluid order in terms 
of the return probability of a random walker. Using 
a recent rigorous result for the latter, we can establish 
bounds on the former for an any site of the cluster, pro- 
vided that we are away from threshold on the percolating 
side (p < p c ). 

Consider a random walker placed on the percolation 
cluster, that begins its walk at site i. The differential 
equation describing the time evolution of the probability 
of this walker being found at time t at position j (Pi(j,t)) 
is: 

d t P i {j,t) = -T jk P i {k,t) (18) 

where an appropriate choice of the stepping rate has been 
made. Thus, the probability of the particle returning to 
the site i where it started from is simply related to the 
eigenvalues and eigenfunctions of the Laplacian operator 
on the percolation cluster: 

Pi(i,t) =J2xl(*)eM-Kt) (19) 

n 

The root mean square (RMS) phase fluctuation ampli- 
tude at a site i ijTUI) can be simply expressed in terms of 
this return probability Pu{t) as 

< $ >= r (20) 

JO V7rf 

As a simple check of this formula, consider the case 
of a linear lattice, uniform, nearest-neighbor hopping 
Jij = Ji.i+i = J ■ Then, the walker just performs a 
random walk in one dimension. If one takes the naive 



continuum limit of (|18|) . it is not hard to see that the 
return probability decays like as t — > oo. Apply- 

ing 1)20(1 . this implies the usual logarithmic divergence in 
the RMS phase fluctuations that destroys LRO in one 
dimensional superfluids 17 . 

For some percolation clusters, Remy and Mathieu have 
found an upper bound for the asymptotic decay of the 
return probability.— They considered site-percolation on 
the square lattice and looked at strictly p < p c : the ran- 
dom walker has an infinite connected cluster to explore. 
They showed rigorously that for almost all cluster real- 
izations, as t — > oo, 

Pi(i,t)<l/t,, (21) 

where, i refers to any site of the cluster. In other words, 
the asymptotic decay of the heat kernel on all sites of 
almost all cluster realizations, decays as fast as it does in 
the undiluted two-dimensional lattice. Using (|20() . this 
implies that for any site on the infinite cluster for p <p ci 
the contribution to RMS phase fluctuations from asymp- 
totically long time-scales is indeed finite, and bounded. 

The two main results of this section suggest that long- 
range order has stability to weak quantum fluctuations. 
We established these results by making the spin-wave 
approximation to linearize the equations of motion. In 
the following section we introduce a complementary ap- 
proach (the slow-blob approximation) that treats the 
non-linear dynamics beyond the spin-wave approxima- 
tion. 



III. SLOW BLOB APPROXIMATION TO O(N) 
MODEL AT p — p c 

Within the slow-blob approximation to the backbone, 
we proceed in two stages. In the first stage, we treat blob 
degrees of freedom as parameters rather than dynamical 
variables, and we solve the link problem for the ground- 
state energy E({Li,fii}). Since links connect only to 
blobs (Fig. 2]) , the link problem reduces to that of inde- 
pendent links. In the second stage, we treat the ground- 
state energy E({Li, Hi}) as an adiabatic potential energy 
for the blobs, and solve for the blob ground-state. 

An analogous procedure arises in the study of molec- 
ular vibration spectra, where it is called the Born- 
Oppenheimer approximation^ In that case, the slow 
degrees of freedom, analogous to blobs, are the mas- 
sive ion cores. The fast degrees of freedom, analogous 
to links, are outer shell electrons. Corrections to the 
Born-Oppenheimer approximation are suppressed by the 
ratio m/M ss 10 -3 , of electron to ion mass. 

Here, the relevant control parameter is the ratio u/U of 
blob charging energy to bare charging energy. If the blobs 
are internally ordered then one has the simple relation 

u=-, (22) 
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FIG. 1: The backbone of a percolation cluster is the union of 
self-avoiding walks on the cluster. It consists of blobs (ovals) 
connected by links. The internal structure of a blob (inset) 
shows that blobs consist of bonds which may be removed with- 
out disconnecting the cluster. In contrast, links consist of 
bonds which, if removed, disconnect the cluster. 



TABLE I: Correlation-length v, backbone db, and blob-size 
a exponents for percolation in various dimensions d. For ap- 
proximate results, number in parentheses indicates standard- 
deviation of last digit. 



d 


2 


3 


> 6 a 


V 


4/3 6 


1.12(1)° 


1/2 


dt 


1.6431(6)° 


1.74(2)° 


2 


a 


0.5435(2) 


0.487(10) 






"M ean -field theory works for d > 6 and gives exact exponents, 
Rcf. 23. 
6 Exact. see Ref.l2l 
c Ref.|2g 

d Ref.Ea 



where, n is the number of sites in a given blob. Now, on 
average, blobs are very big: < n >« ]^d b -i/^^ wuere £, 
is the side-length of finite lattice See Table H] for es- 
timates of v and d&. The bottom line is that, in two 
and three dimensions dt — \jv > 0. For dt > 1/z/, 
the blobs attached to the ends of a given link, have, 
on average, a thermodynamically small charging energy: 
u ~ U / L db ^ 1 / U . Further, even the smallest blobs must 
contain at least four sites, so, even in the worst case the 
blob charging energy u/U as .25 is reasonably small. 

This computation treats the blobs as internally or- 
dered, which is a reliable approximation. If the blobs 
were not internally ordered, then the site-averaged RMS 
phase fluctuations < 4> 2 > pcrc would not be finite. In sec- 
tion II, we showed in the spin-wave approximation that 
the fluctuations are finite (See l|17l) ). Thus the blobs 
must be internally ordered, and our computation in the 
preceding paragraph can be trusted. 

Following the idea of the slow-boson approximation, 
we consider the link problem with the blobs treated as 
parameters and determine P(j) for j <C J. Consider a 
given disorder realization of the otij in the bare Hamilto- 
nian Q. Between two neigboring blobs, the link length 
£ will vary randomly. In fact, the lengths have an expo- 
nential distribution: 20,21 



p(£) ~ exp(-^o), 



(23) 



with, for example, £n w 2.7 for square lattice 
percolation^ If we impose twisted boundary conditions 
with ni ^ fi2, the ground-state energy of the link in- 
creases by an ^-dependent amount: 



E{n\,n2) ~ -j(£)ni ■ m, 



(24) 



where, the "stiffness" j(£), depends on the link length, 
£ and, implicitly, on the parameters of the bare model 
J,U in (QJ. From this we can deduce the distribution 
of blob exchange couplings P(j) appearing in J5J) by a 
simple change of variables: 



dj(t) 



df 



(25) 



where, £ is the link length satisfying j(£) = j. 



It remains to compute the stiffness of a link. Fortu- 
nately, Cardy has computed j(£)m& For N = 2, there 
are two qualitatively different cases to consider: large U 
(U > it 2 J/4) and small U (U < n 2 J/4). For large U, the 
link breaks into segments of length £ with the property 
that rotors in different segments are essentially uncor- 
rected. In this "short-range" phase, the stiffness drops 
off exponentially with link length: j(£) ~ exp(— £/£). In 
contrast, for small U, the stiffness falls off only as £~ x : 
j(£) — * k/£. For N > 3, there is only one case to consider, 
since for all values of U/J, the link is in a short-range 
phase with exponentially small stiffness^ To summarize, 
the stiffness depends on link length as follows: 



k/t, N = 2, (large*/) 



N > 3,andiV = 2, (smaJlfJ). 



(26) 



Plugging this into (|25|l gives the link strength distribution 
P(j) 



P(j) 



e -*/loj j iV = 2, (large*/) 



N> 3,andiV = 2, (small*/). 



(27) 



The distribution R(u) follows directly from the dis- 
tribution of blob sizes n. For large n, the fraction of 
blobs with size n, scales like n _2+Q , with < a < 1, 
given in terms of backbone dimension, db and correlation 
length exponent v via a = 1 — \/{dbv)^ See Table H] 
for estimates of a. Thus, we have determined the charg- 
ing energy distribution entering , using 122|) to change 
variables from blob size n to charging energy u: 



R(u) oc u 



(28) 



where, u <ti U is the charging energy of a given blob. We 
have now completely specified the blob problem J3J. 

The goal of this model J2J is to understand the compe- 
tition between two possibilities: cither very large blobs 
tend to "anchor" order across one-dimensional links, as 
has been suggested to explain ordering in the diluted 
2D Heisenberg model, or else fluctuations within one- 
dimensional links destroy LRO. Both possibilities are re- 
alized in our model for different values of the initial dis- 



6 



tributions. For parameters chosen to reflect actual per- 
colation clusters in 2D, our renormalization group calcu- 
lation, described in section IV, finds LRO. 

The next step is to see how this power-law distribution 
of blob sizes enhances ordering and leads to new fixed 
points in a real-space renormalization group calculation 
for disordered interacting bosons in ID. 



IV. RENORMALIZATION GROUP FLOW OF 
0(2) MODEL FOR LARGE U 

The idea of the real-space renormalization group for 
the 0{2) rotor in ID is to successively integrate out either 
the largest Josephson coupling or the largest charging 
energy (See Appendix A for a description of the elemen- 
tary renormalization group step and its generalization 
to the 0(N) case). This iterative procedure generates 
flow equations for the distributions of charging energy 
and Josephson coupling. The capacitance distribution is 
given as a function / (£) of the scaled variable ( = Cl/u— 1, 
where f2 is an upper cutoff of energy: fl — maxijiti, ji}, 
and the Josephson coupling distribution is given as a 
function g([3) of the scaled variable (3 — log(£l/j). 

The analysis of one-dimensional models begins with 
the RSRG flow equations^ as a function of energy scale 
r for the charging energy distribution /(£, T) and Joseph- 
son coupling distribution g(f3, T): 



dj_ 



dg_ 
dT 



pOC pOO 

= 9o / /(Ci)/(C 2 )5(Ci + C 2 + i - dCi d( 2 

Jo Jo 

+ (l + 0§£ + (/o + l-So)/. 

p oo poo 

= f / dfo dp 2 g((3 1 )g(p 2 )S(f3 1 + fa - [3) 
Jo Jo 

+ || +5(9o-/o). (29) 



Here T = log(Sl/ /f2) tracks the progress of the renor- 
malization flow (fli is the largest coupling in (J2J) before 
renormalization), fo = /(O) and go = g(0). We will even- 
tually discuss the behavior of the coupled / and g equa- 
tions, but for now, consider go as a constant parameter 
in the / flow equation. 

It is possible to get immediate insight into the phase 
diagram found for exponential distributions^ by aver- 
aging both sides of the flow equation. The average 
(£) = C/(C) d(, if it exists as for the exponential dis- 
tributions considered in by Altamn et. al>i&, evolves as 



dT 



<C)=.9o(i + (C)) + /o(C)-i-(C). 



(30) 



In this form, we can understand the appearance of a tran- 
sition at go = 1 for the class of distributions studied 
in |l6j : for that one-parameter family of distributions 



fa(C) 



ae 



-to, 



, so that fo = a, we have (£) = 1/a so 



The flow equation for the mean capacitance, if the func- 
tion / is initially of the form / a , is 



(C) - go(l + 1/a) - 1/a = .90 + (<?o - I) /a. (32) 



The replacement <5(Ci + C2 + I — £) — * <5(Ci +C2 — C) made 
inl£ corresponds to replacing g (l + {Q) by go(0 in i|3U|l . 
With this replacement, an / distribution initially in the 
one-parameter family remains in this one-parameter fam- 
ily even for nonzero go, and the average charging energy, 
a, flows as follows^ 



da 

^ = (!-*)«. 



(33) 



This equation is written in terms of a rather than fo in 
order to stress that it is specific to exponential distribu- 
tions. It is now shown that power-law distributions have 
a different RG equation with different qualitative behav- 
ior, leading to a significantly different phase diagram. 

There are exact fixed points of the RG flow equation 
that are not contained in the family of exponential distri- 
butions that has been studied previously. Power-law dis- 
tributions are needed both to describe these fixed points 
and to model the fractal distribution of blob sizes on the 
percolation cluster [See Eq. [28]. Setting go = 0, we ob- 
tain a one-parameter family of static solutions of l|29|) 



M0 = 



(i + C) 1+fc 



(34) 



These are well-behaved decreasing normalized distribu- 
tions for all b > 0. For b > 1, the mean (() converges 
and satisfies 



/o(C) - <C) - 1 = 



(35) 



so that, as required, the mean is constant in time. Hence 
there is a one-parameter family of fixed points that lies 
outside the one-parameter family f a considered in [lfij . 
which contains no fixed points at finite a. Now if the 
function / is originally of the form the mean (£) is 
increasing if go > 0, suggesting that the system can flow 
to large capacitance (small charging energy) if go > 0. 

More explicitly, we simply evaluate the right-hand side 
of the flow equation 129|) using our power-law distribu- 
tions /b(C), and find that the value of the distribution 
at large capacitance increases. That is, the distribution 
flows to large capacitance assuming that it starts out with 
a power-law form. For a power-law form, the integral 
appearing in l|29|l has the following simple, asymptotic 
behavior: 
1 

/6(Ci)/&(C 2 )5(Ci + C2 + 1 - C) dCi dC 2 -> 2/6(0, 

(36) 

for C ^ 1- Using this fact, and the fact that at go = 0, 
power-laws are fixed points of the flow, one finds that the 
fb(C) distributions flow as follows at large capacitance: 



/o(C)-l-<C> =-!/<* = -V/o. 



(31) 



dfb 
dT 



9oMQ. 



(37) 
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For go > 0, the right-hand side is positive, and, as the 
result, the value of the distribution at large values of ( 
increases under the flow. Thus, we see explicitly the ten- 
dency for power-law distributions to flow towards large 
capacitance (small charging energy). 

So far we have seen that the mean (£), assuming 
go = 0, is constant for the family fb and decreasing 
for the family f a . There are also well-behaved distri- 
butions for which /o(C) — (C) ~ 1 > 0: an example is 
/(C) = 10(1 + VC)~ 6 , for which f = 10 and (0 = 1. 
However, one might wonder if the initial increase or de- 
crease of the mean (£) is a useful guide to the eventual 
behavior of the RG flow. To answer this question, Ap- 
pendix B constructs additional soliton solutions, which 
include the fixed points fb(C) as a limiting case, and for 
which the go = flow equation is exactly solvable. One 
family of solutions are related by a logarithmic rescaling 
of variables to solutions discussed by Fisher in the con- 
text of spin chainsi^ The construction exhibits an inter- 
esting cumulant property of the nonlinear flow equation 
at go = 0. 

The previous paragraphs constructed a new set of 
fixed-point solutions, the /&(C) in 1|34|) . to the functional 
RG equation with go = 0. Appendix B gives other solu- 
tions to the go = equation for /, but no full solution to 
the coupled RG equations with fo and go both nonzero 
has been found. Hence some other means of projecting 
the infinite-dimensional renormalization group flows to 
a finite-dimensional subspace is necessary. This section 
gives a prescription for this pro jection that both justifies 
the phase diagram found in [T|| for exponential distribu- 
tions and gives a new phase diagram once power-law tails 
in capacitance are allowed. The RG flows superficially 
depend on a cutoff introduced in the projection process, 
but the phase diagram and location of the critical point 
are cutoff-independent. 

Let < a < 1 be some arbitrary cutoff, and define 
< W a < 1 as the integral 



W a = 



/(C) d(. 



(38) 



The evolution equation for W a is 
dW a 



dl 



= foW a -(l + a)f(a)+g (l-W a ). (39) 



Note that for g = 0, the fixed-point solutions fb(() have 
W a = 1/(1 + a) b so that dW a /dT = correctly. In a mo- 
ment we will consider the flow generated by this equation 
within the fb, but first consider the exponential solutions 
/(C) = /ocxp(— /oC)- Substituting into l|3*9")l gives the 
projected flow 



dfo 



fo ~ 3o- 



3 a/o _ ^ 



dT a 
When fo is small, this equation becomes 

dfo 



dr 



fo — fo9o, 



(40) 



(41) 



which is exactly the RG flow obtained in 0] under an ap- 
proximation valid in this small- fo limit. The same equa- 
tion is also obtained taking a — > for any fo. For small 
fo, the superfluid transition in these equations at go = 1 
has been discussed thoroughly by Altaian et alii 6 - The 
critical point fo — 0, go — 1 corresponds to the transition 
between "stiff" and "floppy" regimes of the classical XY 
model. 26 This RG flow is shown in Fig. 1. 

Now consider instead the flow within the power-law 
solutions f . For g ^ 0, the evolution equation (f3~9"|l de- 
fines a projected flow within the one-parameter subspace 
fb- since f = /(0) = b, we have 



dfo 
dT 



-,9o 



l-W a 



W a log(l + a) 



= -.90- 



exp(/ log(l + a)) - 1 



log(l + a) 



(42) 

Again taking the limit of small fo and any a, or small a 
and any fo, we have 



dfo 
dT 



-fo9o- 



(43) 



Note that the right-hand-side is the same as in the pre- 
viously obtained flow equation for coupling strength 16 



dgo 
dT 



= ~fo9o 



(44) 



so that fo — go is constant and the flows are as shown in 
Fig. 2. This coupling strength flow equation describes 
the flow within a power-law family of distribution of cou- 
pling strength P(j) oc j 90 ^ 1 . As argued in the previous 
section, such a distribution enters (0 for the 0(2) rotor 
model with large U [See Eci. l27|) ]. In the next section, we 
discuss the small U case. 

Note that now there is no phase transition at the value 
go = 1; instead there are two lines of fixed points, one 
with fo = and one with go = 0. The entire classical line 
fo = is stable, including both the "stiff" and "floppy" 
regimes of the classical XY model. Now, we computed 
the initial conditions for this flow (See (|27|l and (|28|) ~): 



/o(r = 0) = l-a, 

<? (r = o) = 



(45) 



where, £ is the correlation length of a 0(2) rotor chain 
with large charging energy U, £o is the average length 
of links of red sites, and, a = 1 — l/(dbv) is the blob 
exponent tabulated in Table [I] Since fo — go is a renor- 
malization group invariant of the flow in (|43I44|) . we can 
determine the fixed point towards which the flow tends 
in a simple way. Namely, the flow tends towards the 
ordered, classical line fo = provided that initially we 
have go(T = 0) > /o(r = 0). That is, to have long-range 
order, we must tune U / J close enough to the Kosterlitz- 
Thouless transition of a one- dimensional, clean 0(2) ro- 
tor model at U/J ~ tt 2 /4, where £ diverges, such that 
the following criterion is satisfied: 



S>to(l-a) 



(46) 
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fo 




FIG. 2: Renormalization-group flows projected to two- 
parameter plane, after |T3| . for capacitance distributions with 
exponential damping at large capacitance. 



fo 




go 



FIG. 3: Renormalization-group flows projected to two- 
parameter plane, for capacitance distributions with power- 
law damping at large capacitance. Note the new line of fixed 
points at <jo = and the absence of a phase transition at 
(9d,/o) = (l,0). 



Further, when this criterion is not satisfied, the flow 
tends towards the insulating line of fixed points go = 0. 
Thus, the blob hamiltonian J2J admits both insulating 
and superfiuid ground-states. When U/J is small enough 
to satisfy criterion ()4fci[l. the renormalization group flow 
tends to the superfiuid ground-state. The possibility of 
an ordered phase stable to nonzero charging energy is 
novel, and supports the existence of order for the 0(2) 
rotor model on percolation clusters via the connection 
discussed in the Introduction. 



V. 0(2) ROTOR FOR SMALL U 

In this section, we consider the special case of the 0(2) 
rotor with small U, at percolation threshold, that is, |QJ 
with N = 2. We will argue that this case is only slightly 
different than the large U case considered in the pre- 
vious section. Namely, we will show using a spin-wave 
argument and a real-space renormalization group argu- 
ment, that this case exhibits long-range order regardless 
of the value U/J, so long as U is small. Recall that, by 
contrast, for large U, one must satisfy the criterion l|4t)|) 
for the ground-state to exhibit long-range order. 

As in the previous section, we begin with the slow-blob 
approximation J5J) . In this approximation, the novelty of 
the 0(2) rotor with small U arises from the "activated" 
form of the exchange distribution P(j) in (1271) : 



P{j) 



,-k/ioj 



(47) 



In other words, the power-law charging energy distribu- 
tion H28fl holds for this case too. Only the exchange distri- 
bution changes compared to the large U case. Notice that 
the new, activated exchange distribution is significantly 
less broad than the power-law distributions considered in 
the previous section. 

Let us, to first approximation, neglect completely the 
spread in the exchange distribution, for this activated 
form. This is a drastic approximation which we will im- 
prove upon below by considering the flow of the distri- 
bution under real-space renormalization. However, the 
activated form (|47|l is quite narrow, and treating it as 
completely concentrated at its mean-value will turn out 
to be a consistent approximation. In this approximation, 
we may use the following exact result for the asymptotic 
decay of the (appropriately defined) return probability of 
a random walker: 



( Pi (i,t))^(—)T=*7* 



(48) 



Alexander et. ali^I show that this result obtains for uni- 
form exchange distribution and power-law charging en- 
ergy distribution R(u) cx u~ a , with < a < 1. Recall, 
that this range of a obtains for percolation in two and 
three dimensions (See table P). Using l|2U|) . which applies 
to , the asymptotic result in 148(1 implies a finite RMS 
phase fluctuation of the blob phases fa, where, the blob's 
rotor variable in @ is given by n$ = (cos(fa) , sia(fa)) , 
for N = 2 rotor. Thus, in the approximation of uniform 
exchange energy, spin-wave theory predicts that © ex- 
hibits an ordered phase stabilized by the small-£7 part of 
the charging energy distribution. 

Let us now consider the actual distribution (|47() and 
see how it flows under real-space renormalization. As 
in the previous section, we gain immediate insight into 
the flow by averaging both sides of the flow equation 
(j5SJ. The average (/?) = J °° f3g((3) df3 exists for activated 
distributions l|47|l and evolves as follows: 



clT 



(/3> = -l + (<?o + /o)</3>. 



(49) 
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Recall, that j3 = log(f2/j) is the scaled Josephson ex- 
change and g((3) = P(j — Qe~P)Qe~P is the distribution 
of /?, where, P(j) is given by 147(1 . More directly, the 
activated ansatz (|47|) implies a distribution of j3 of the 
following form: 



;(/3,r)= 5o (r)e^( r )e a -/ 3+9 o(r )j 



(50) 



Here, T = log(il//57) is the flow parameter, where, 
ill = m a x {ji,Ui} is the largest coupling in @ before 
any renormalization, and f2 < 0/ is the largest coupling 
at the current stage of renormalization. The parameter, 
5o(r) initially satisfies go (T = 0) = k/£o&i, ■ At later 
times, we can evaluate how 30 (r) flows by first computing 
the expectation appearing in l|49|) : 



(/3}=e 9o(r) £i(So(r)). 



(51) 



Here, the exponential integral appears: E\{go) = 
J°°d7e-T/7. Substituting |J5TJ into |@|J gives the pro- 
jected flow of go(r) : 



dgo = -l + (fo+go)e 9o E 1 ( g ) 
dr -l+g e3oE 1 (g ) 



-go- 



(52) 



When go is small, the flow simplifies to the following 
dgo 



dT 



go + /o.9o In 3o 



(53) 



where, the omitted terms are higher order in go. The 
right-hand side of (|53[1 is initially positive, so go grows un- 
der renormalization. Recall, that go is the weight of the 
distribution at the large values of the exchange strength 
j. The projected flow equation implies that our ac- 
tivated form H47JI flows towards a narrower distribution 
concentrated at the largest values of j. This justifies 
the spin-wave treatment at the beginning of this section 
(Eq. Q48JI and paragraph containing it), which assumed a 
narrow distribution of j. 

Thus, the conclusion suggested by spin- wave theory 
holds even when we allow the exchange distribution to 
start with an activated form. Since this is precisely 
the form of distribution satisfied by the 0(2) rotor with 
U/J < 7r 2 /4, the model must have finite RMS fluctu- 
ations of the blob rotor variables Hi. In the slow-blob 
approximation, this implies that the percolation cluster 
as a whole exhibits long-range order for the 0(2) rotor 
model with U/J < 7r 2 /4. 



VI. CONCLUSION 

To conclude, the question of whether the O(N) (N > 2) 
quantum rotor model can have an ordered phase on the 
percolation cluster was studied via a spin wave calcula- 
tion and via a real-space renormalization group calcula- 
tion. The result for these models may be summarized 
in the following way which also tracks the details of the 



calculation. We associate with the incipient percolation 
cluster an effective spatial dimensionality in order to dis- 
cuss the question of ordering. The correct spatial dimen- 
sionality for this purpose is the fracton dimension d s , 
defined in equation l|16|) from the spectral density of the 
Laplacian on the percolation cluster. This dimensionality 
d s is distinct from the fractal dimension of the percola- 
tion cluster, and happens to be numerically very close to 
4/3 for percolation in any dimension. The possibility of 
long range order can now be determined by simply com- 
paring this dimension with the lower critical dimension 
arising from the Hohenberg-Mermin- Wagner theorem. 17 
Thus, for the zero temperature quantum models where 
d s + 1 > 2 we have the possibility of long range order 
for percolation in any dimension, while at finite temper- 
atures, these models are always disordered since d s < 2. 

The presence of weak links compelled us to investi- 
gate further the possibility of long range order. We "in- 
tegrated" out the weak links, generating a low-energy 
description [See Eq. J2J], amenable to real-space renor- 
malization. The renormalization group flow [See Fig. [3] 
and Eq. (|43|) . l|44|) ] demonstrated a competition between 
two natural length-scales associated with the links [See 
Eq. 14t)|) ]. On the one hand, the average length of a link 
£o, and on the other, the correlation length within a long 
link £. The final result may be stated simply: When 
U/J is small enough that £ > ^o(l — then correla- 
tions "jump" across the weak links and create long-range 
order. The renormalization-group flow supporting this 
conclusion depends crucially on the fractal structure of 
the percolation backbone. 

These considerations are summarized in the phase di- 
agrams in Fig. 0] In Fig. 0Ja) , the existence of a finite 
superfluid order parameter at the percolation transition 
for small U/J leads to the vertical phase boundary. In 
Fig. Efb) , we simply point out that the well-known ab- 
sence of a Mott insulating state for incommensurate fill- 
ing at zero dilution persists under dilution until percola- 
tion threshold. Finally, the absence of superfluid order 
at finite temperatures on percolation clusters in any di- 
mension D > 2 gives rise to the phase diagram shown in 
Fig. He). 
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FIG. 4: The phase diagram for bosons with onsite repulsion 
on a diluted lattice as a function of the dilution fraction p. 
The percolation transition occurs at p c . (a) Integer filling 
at zero temperature: The Mott insulating phase arises when 
charging energy U is large compared to Josephson coupling J. 
(b) Away from integer filling, the Mott phase is not present 
for any U /J. (c) Finite temperature transition for bosons out 
of the superfluid phase. 



APPENDIX A: RG OF RANDOM O(N) CHAIN 

Consider the O(N) rotor chain with random exchange 
and "charging" given in Eq. ■ We have two possible rg 
steps. First, if the largest energy is an exchange Ji, then 
we lock the two rotors together, and estimate the effective 
charging energy of the combined rotor. The result is that 
inverse charging energies add: 



u- 



u; 



(Al) 



Second, the largest energy may be a charging energy Q = 
Ui. In this case, we put this rotor in the rotationally 
invariant ground-state, and treat the neighboring rotors 
as a classical magnetic field h = Jifii + Ji+\fli + 2- For 
small the ground-state energy is, to lowest order, 

2h 2 /N(N — 1)0. Expanding this, yields a new effective 
coupling between the neighboring rotors: 



J 



4J ( J 4 



n(n- i)fi 



(A2) 



The first step is identical to 0(2) case. The second step 
is virtually identical. If we measure exchanges on a log- 
arithmic scale, Pi = log(S!/Ji), then the O = U step 
becomes 



/3 = A+/3 i+1 +l g(iV(iV-l)/4). 



(A3) 



As always, in asymptotic analysis of these kind of flow 
equations, we can neglect additive constants in the RG- 
step, since we are looking for broadly distributed fixed 
points^ In conclusion, we have obtained the elementary 
RG-step for the O(N) rotor chain with strong disorder. 



APPENDIX B: ADDITIONAL SOLUTIONS OF 
THE RSRG FLOW EQUATION 

The goal of this section is to find Independent exact 
solutions of the capacitance flow equation i|29|) that ver- 
ify the statements in Section III derived using projected 
RG flows. Although we have not been able to solve the 
nonlinear PDEs for arbitrary initial data, the exact solu- 
tions now discussed are examples supporting the qualita- 
tive behavior found using a moment equation in Section 
III. Let us look for a solution in terms of a distribution 
0(6, r) over the one-parameter family fb(C) : 



/(C,r) 



<t>(b,T)f b (0db= / 0(6, r) 



1 



db. 



1+b 

(Bl) 

Now /o(r) = /(0,T) = f °° 6^(6, T) db, i.e., the average 
(6)0 taken over the distribution <j)(b, T). Our goal will be 
to find solutions of the evolution equation for 0(6, T): 



d</>(b,T) 
dT 



-b + 



0(6 2 )6 2 d6 2 U(6, T). (B2) 



Now assume that 0(6, T) is nonnegative (it is then auto- 
matically normalized to 1) so that it can be regarded as 
a probability distribution. Its normalization is constant, 
and its mean (equal to /o) evolves via 



d(b) 4 

dr 



(B3) 



In terms of the mean /zi and the variance /i 2 , this equa- 
tion is compactly written = — /i 2 . A bit of algebra 
confirms that a similar evolution holds for all cumulants: 



dT 



(B4) 



This provides a way to find soliton-like solutions. One 
class of solutions starts from the generalized Poisson dis- 
tribution with two parameters A and xq, 



0(6) = ]T<S(6-fc C o) 



X k e- 



(B5) 



which has cumulants [i n = Axo™. This distribution has a 
finite probability of 6 — 0, which corresponds to exactly 
charging energy; below we find a different soliton solution 
without this property. The discrete distribution l|B5|l ap- 
proaches a continuous Gaussian centered at 6 = 60 if we 
take the limits xq — > and A — > 00 with Xxq — 60 . Then 
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the cumulant evolution equations are all solved if xq is 
constant and A(T) solves 

Y = -x X => A(r) = \(0)e- x ° r . (B6) 

In the limit of small xo with xqA(O) fixed, we recover the 
fixed-point solutions found previously. 

A second family of solutions starts from the gamma 
distribution with parameters a and 9: 



0(6) = 



b a-i e -b/e 



oT 



if 



(B7) 



T{a)6 a 

The cumulants of the above distribution are 

K n =T(n)a6 n . (B8) 
So, if the only T dependence is in 9, 



T(n + I)a6 n - L ^ = -a9 n+1 T(n + 1) = -K n+1 

(B9) 

(BIO) 



It is simple to confirm explicitly that the resulting form 
for /, 



/(C,r) 



a6(T) 



(i + C)(i + 0(r)io g (i + O) 



a+1 ' 



(Bll) 



is normalized for a > and satisfies the evolution equa- 
tion. For this solution, / (r) = a/(T + 9(0)- 1 ), indicat- 
ing that this charging energy distribution is marginally 
irrelevant even as go — > 0, well away from the classical 
transition. 

We also note that the gamma solutions have normal- 
ized variance K2/(«i) 2 = l/c*, which is a constant of the 
motion. As a — > 0, the distributions become extremely 
broad, and the RG is expected to be valid in this limit. 

These two families of solutions suggest that having a 
finite probability of zero charging energy, as in the Pois- 
son case, leads to an exponentially rapid flow of the entire 
distribution to zero charging energy: the order is strongly 
"nucleated" by the classical rotors with zero charging en- 
ergy. The gamma solutions also flow to zero charging 
energy, proving that such flow is possible for regular dis- 
tributions, but this solution flows only as a power-law. 
Note that the flow equation for 9 is exactly that for a 
marginally irrelevant operator. 
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